(N 

o 

O 

Q 

(N 



O 

in 



O 

o 



> 
o 

in 

(N 



X 



Janus fluid with fixed patch orientations: theory and simulations 
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We study thermophysical properties of a Janus fluid with constrained orientations, using analytical techniques 
and numerical simulations. The Janus character is modeled by means of a Kern-Frenkel potential where each 
sphere has one hemisphere of square- well and the other of hard-sphere character. The orientational constraint 
is enforced by assuming that each hemisphere can only point either North or South with equal probability. 
The analytical approach hinges on a mapping of the above Janus fluid onto a binary mixture interacting 
via a "quasi" isotropic potential. The anisotropic nature of the original Kcrn-Frcnkel potential is reflected 
by the asymmetry in the interactions occurring between the unlike components of the mixture. A rational- 
function approximation extending the corresponding symmetric case is obtained in the sticky limit, where 
the square-well becomes infinitely narrow and deep, and allows a fully analytical approach. Notwithstanding 
the rather drastic approximations in the analytical theory, this is shown to provide a rather precise estimate 
of the structural and thcrmodynamical properties of the original Janus fluid. 



I. INTRODUCTION 

Janus fluids refer to colloidal suspensions formed by 
nearly spherical particles with two different philicitics 
evenly distributed in the two hemispheres.—^ Under typ- 
ical experimental conditions in a water environment, one 
of the two hemispheres is hydrophobic, while the other 
is charged, so that different particles tend to repel each 
other, hence forming isolated monomers. On the other 
hand, if repulsive forces are screened by the addition of 
a suitable salt, then clusters tend to form driven by hy- 
drophobic interactions^ 

This self-assembly mechanism has recently attracted 
increasing attention due to the unprecedented improve- 
ment in the chemical synthesis and functionalization of 
such colloidal particles, that allows a precise and reliable 
control on the aggregation process that was not possible 
until a few years ago^ From a technological point of view, 
this is very attractive as it paves the way to a bottom- 
up design and engineering of nanomaterials alternative 
to conventional top-down techniques^ 

One popular choice of model describing the typical du- 
ality characteristic of the Janus fluid is the Kern-Frenkel 
models This model considers a fluid of rigid spheres hav- 
ing their surfaces partitioned into two hemispheres. One 
of them has a square- well (SW) character, i.e., it attracts 
other similar hemispheres through a SW interaction, 
thus mimicking the short-range hydrophobic interactions 
occurring in real Janus fluids. The other part of the sur- 
face is assumed to have hard-sphere (HS) interactions 
with all other hemispheres, i.e., with both like HS as well 
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as SW hemispheres. The HS hemisphere hence models 
the charged part in the limit of highly screened interac- 
tions that is required to have aggregation of the clusters. 

Although in the present paper only an even distribu- 
tions between SW and HS surface distribution will be 
considered (Janus limit), other choices of the coverage, 
that is the fraction of SW surface with respect to the 
total one, have been studied within the Kern-Frenkel 
models In fact, one of the most attractive features of the 
general model stems from the fact that it smoothly in- 
terpolates between an isotropic HS fluid (zero coverage) 
and an equally isotropic SW fluid (full coverage).^ 

The thermophysical and structural properties of the 
Janus fluid have been recently investigated within the 
framework of the Kern Frenkcl model using numerical 
simulations^^ thus rationalizing the cluster formation 
mechanism characteristic of the experiments £ The fluid- 
fluid transition was found to display an unconventional 
and particularly interesting phase diagram, with a re- 
entrant transition associated with the formation of a clus- 
ter phase at low temperatures and densities While 
numerical evidence of this transition is quite convincing, 
a minimal theory including all necessary ingredients for 
onset of this anomalous behavior is still missing. Two 
previous attempts are however noteworthy. Reinhardt 
et aZAi introduced a van der Waals theory for a suitable 
mixture of clusters and monomers that accounts for a re- 
entrant phase diagram, whereas Fantoni et alr^^ devel- 
oped a cluster theory explaining the appearance of some 
"magic numbers" in the cluster formation. This notwith- 
standing, the challenge of an analytical theory fully de- 
scribing the anomaly occurring in the phase diagram of 
the Janus fluid still remains. 

The aim of the present paper is to attempt a new route 
in this direction. We will do this by considering a Janus 
fluid within the Kern-Frenkel model, where the orienta- 
tions of the SW hemispheres are constrained to be along 
either North or South, in a spirit akin to Zwanzig model 
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for the isotropic-nematic transition in liquid crystals^ 

Upon observing that under those conditions, one ends 
up with only four possible different interaction (North- 
North, North-South, South- North, and South-South), 
this constrained model will be further mapped onto a 
binary mixture interacting via a "quasi" isotropic poten- 
tial. Here the term "quasi" refers to the fact that a cer- 
tain memory of the original anisotropic Kern-Frenkel po- 
tential is left: after the mapping, one has to discriminate 
whether a particle with patch pointing North ( "spin- up" ) 
is lying above or below that with a patch pointing South 
("spin-down"). This will introduce an asymmetry in the 
unlike components of the binary mixture, as explained 
in detail below. In order to make the problem tractable 
from the analytical point of view, the particular limit of 
an infinitely narrow and deep square- well (sticky limit) 
will be considered. This limit was originally devised 
by Baxter and constitutes the celebrated one-component 
sticky-hard-sphere (SHS) or adhesive Baxter model.— By 
construction, our model reduces to it in the limit of fully 
isotropic attractive interactions. The latter model was 
studied within the Percus-Yevick (PY) closure 1 ^ in the 
original Baxter work and in a subsequent work by Watts 
et alrH The extension of this model to a binary mixture 
was studied by several authors^ - — The SHS model 
with Kcrn-Frenkel potential was also studied in Ref. [23l . 
via a virial expansion at low densities. 

A methodology alternative to the one used in the 
above studies hinges on the so-called "rational-function 
approximation" (RFA) ; 24 i 25 and is known to be equiv- 
alent to the PY approximation for the one-component 
SHS Baxter models and for its extension to symmetric 
SHS mixtureSi 18 i 22 ' 24 The advantage of this approach is 
in the fact that it can be readily extended to more general 
cases, and this is the reason why it will be employed in 
the present analysis to consider the case of asymmetric 
interactions. We will show that this approach provides a 
rather precise estimate of the Janus fluids with up-down 
orientations by explicitly testing it against Monte Carlo 
(MC) simulations of the same Janus fluid. 

The remaining part of the paper is envisaged as fol- 
lows. Section [XT] describes our Janus model and its map- 
ping onto a binary mixture with asymmetric interactions. 
It is shown in Sec. IIIII that the thermophysical quantities 
do not require the knowledge of the full (anisotropic) pair 
correlation functions but only of the functions averaged 
over all possible North or South orientations. Section ITVl 
is devoted to the sticky- limit version of the model, i.e., 
the limit in which the SW hemisphere has a vanishing 
well width but an infinite depth leading to a constant 
value of the Baxter parameter r. The exact cavity func- 
tions to first order in density (and hence exact up to 
second and third virial coefficients) in the sticky limit 
are worked out in Appendix [A] Up to that point all the 
equations are formally exact in the context of the model. 
Then, in Sec.|V]we present our approximate RFA theory, 
which hinges on a heuristic extension from the PY solu- 
tion for mixtures with symmetric SHS interactions to the 



realm of asymmetric SHS interactions. Some technical 
aspects are relegated to Appendices IB1 and \C\ The predic- 
tion of the resulting analytical theory are compared with 
MC simulations in Sec. IVI1 where a semi-quantitative 
agreement is found. Finally, the paper is closed with 
conclusions and an outlook in Sec. IVIII 



II. MAPPING THE KERN-FRENKEL POTENTIAL 
ONTO A BINARY MIXTURE 

A. The Kern-Frenkel potential for a Janus fluid 

Consider a fluid of spheres with identical diameters 
a where the surface of each sphere is divided into two 
parts. The first hemisphere (the green one in the color 
code given in Fig.[T|) has a SW character, thus attracting 
another identical hemisphere via a SW potential of width 
(A — I) a and depth e. The second hemisphere (the red 
one in the color code of Fig. QJ is instead a HS potential. 
The orientational dependent pair potential between two 
arbitrary particles /x and v (fj,, v = 1, . . . , jV, where N is 
the total number of particles in the fluid) has then the 
form proposed by Kern and Frenkel^ 

$ (r M „, n M , n„) = c/> H s (»v) + </>sw (?>) * (r 



where the first term is the HS contribution 



^hs (r) 



oo, < r < a, 
0, a < r, 



(2-1) 



(2.2) 



and the second term is the orientation-dependent attrac- 
tive part, which can be factorized into an isotropic SW 
tail 



0sw (r) 



— e, a < r < Act, 



0, Act < r, 
multiplied by an angular dependent factor 



(2.3) 



, _ J 1, if V > and • < 
I 0, otherwise. 

(2.4) 

Here, r M „ = T^/r^, where r M „ = r v — r M , is the unit vec- 
tor pointing (by convention) from particle [i to particle v 
and the unit vectors n M and n„ are "spin" vectors associ- 
ated with the orientation of the attractive hemispheres of 
particles fi and v, respectively (see Fig.Q}. An attractive 
interaction then exists only between the two SW portions 
of the surface sphere, provided that the two particles are 
within the range of the SW potential. 



B. Asymmetric binary mixture 

We now consider the particular case where the only 
possible orientations of particles are with attractive caps 



r 



FIG. 1. The Kern-Frenkel potential for Janus fluids. 



pointing only cither North or South with equal proba 
ity, as obtained by Fig. [T]in the limit n M = z, = — z, 
and with z pointing North. 

Under these conditions, one then notes that the Kern- 
Frenkel potential (|2.1[) - (|2.4[) can be simplified by associ- 
ating a spin i = 1 (up) to particles with SW hemispheres 
pointing in the North z direction and a spin j = 2 (down) 
to particles with SW hemispheres pointing in the South 
— z direction, so one is left with only four possible con- 
figurations depending on whether particles of type 1 lie 
above or below particles of type 2, as illustrated in Fig. [21 
The relationship between the genuine Janus model (see 
Fig. [1]) and the up-down model (see Fig. [3]) is reminis- 
cent to the relationship between the Heisenberg and the 
Ising model of ferromagnetism. From that point of view, 
our model can be seen as an Ising-likc version of the orig- 
inal Janus model. A similar spirit was also adopted in 
the Zwanzig model for the isotropic-nematic transition in 
liquid crystals^ 

The advantage of this mapping is that one can dis- 
regard the original anisotropic Janus-like nature of the 
interactions and re-cast the problem in the form of a bi- 
nary (equimolar) mixture of particles of type i,j = 1,2 
interacting via an asymmetric potential of the form 

= rW HimsW <.; : + <Pji(nj)9(-COs6ij), 

(2.5) 

where cos 6ij = r,j • z (recall our convention r,j = r 3 — r^) 
and 



<Pij(r) = 0hsM 




SW 



(r), if i = 1 and j = 2, 
otherwise. 



(2.6) 



In Eq. ([231) 8(x) = 1 and for x > and x < 0, 
respectively. 

It is important to remark that, in general, fi2(r) ^= 
<P2i(r), as evident from Eq. (|2.6p . Thus, the binary 
mixture is not necessarily symmetric [unless e = or A = 
1 in Eq. (]2.3[) ]. unlike standard binary mixtures where 
this symmetry condition is ensured by construction. In 
the potential (|2.5|) . there however is still a "memory" 






T 



FIG. 2. (Top-left) A particle of type 1 is "below" another 
particle of type 1 providing SW/HS=HS interactions. (Top- 
right) A particle of type 1 is "below" a particle of type 2 lead- 
ing to SW/SW=SW interactions. (Bottom-left) A particle of 
type 2 is "below" a particle of type 1 yielding HS/HS=HS 
interactions. (Bottom-right) A particle of type 2 is "below" 
another particle of type 2 thus leading again to HS/SW=HS 
interactions. 



of the original anisotropy since the potential energy of a 
pair of particles of species i and j separated a distance r 
depends on whether particle j is "above" (cos#,j > 0) or 
"below" (cos 8ij < 0) particle i. In this sense, the binary 
mixture obtained in this way is "quasi" , and not "fully" , 
spherically symmetric. 



Another important point to be stressed is that, while 
the sign of cos 9ij represents the only source of anisotropy 
in the above potential 0y (r^), this is not the case for the 
corresponding correlation functions, which will explicitly 
depend upon the relative orientation cos 9ij and not only 
upon its sign. This applies, for instance, to the pair corre- 
lation functions gij(Y) = 9ij{r',0), as shown in Appendix 
lAl to first order in density in the sticky limit (see Sec. 
IIVj) . As an illustration, Fig. [3] shows the first-order pair 



correlation functions g[\ (r) and gj^ (r) as functions of 
the radial distance r for several orientations 0. 



As our aim is to remove the oricntational dependence 
in the original potential altogether, a further simplifica- 
tion is required to reduce the problem to a simple binary 
mixture having asymmetric correlation functions depen- 
dent only on distances and not on orientations of spheres. 
This will be the orientational average discussed in the 
next section. 
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FIG. 3. (Top) Plot of g[\\r) as a function of r for (a) 9 = 
and 7T, (b) = f and and (c) 6» = f . (Bottom) Plot of 

the regular part of <?i2^( r ) as a function of r for (a) 9 = 0, (b) 
= \ , (c) = f , (d) 9 = and (e) 9 = n. The interaction 
potential is given by Eq. (|2 . 6|) (sketched in Fig. [2}, except 
that the sticky limit with Baxter's temperature r = 0.1 has 
been taken (see Sec. IIV[) . 



III. ORIENTATION AL AVERAGE AND 
THERMODYNAMICS 

A. Orientational average 

Most of the content of this section applies to a mix- 
ture (with any number of components) characterized 
by any anisotropic potential fair^j) = 4>u{— ry) ex- 
hibiting the quasi- isotropic form (|2.5[) . where, in general 
fij(r) ^ V'ji( r ) if * 7^ J- I 11 that case, we note that the 
thermodynamic quantities will generally involve integrals 
of the general form 



dr ff y(r)Jy(r) 



(3-1) 



with 



(r) = F l3 (r)G (cosfl) + i^(r)6 (- cosfl) , (3.2) 

where in general Fij(r) ^ Fji{r) if i ^ j- This strongly 
suggests that one can define the two orientational aver- 
ages g+Ar) and g~Ar) as 



9%{r) =g lJ (r) 



d(cos9)gij(r), (3.3) 



Note that g^ir) = g'jAr), and this suggests the use of 
the notation (r) and g^ (r ) instead of g^ (r ) and g~^ (r) , 
respectively. Taking into account Eqs. (|3.2[) - (|3.4[) . Eq. 
(|3.1[) becomes 

4i = \j dr [9iA r ) F ij W+hitfFji (r)] . (3.5) 
In the case of a double summation over i and j, 

XiXjlij = ^2 XiXj J r/r 77,, ;/•;/•*,, (r) , (3.6) 

where Xi denotes the mole fraction of species i. 



B. Thermodynamics of the mixture: energy, virial, and 
compressibility routes 



We can now particularize the general result p.6[) to 
specific cases. 

In the case of the internal energy J 7 ^ (r) = fa (r) and 
so the energy equation of state can be written aai£ 

u cx = ^p^XiXj J drgij (r) fa (r) 

i,3 

= >■->:, j dvy l0 (r)^{r)e-^\ (3.7) 

i,3 

where u cx is the excess internal energy per particle, p is 
the number density, = 1/ksT (feg and T being the 
Boltzmann constant and the temperature, respectively), 
and Vij(r) = ~ 9 ij{r)e^ ipi A r ) is the orientational average of 

the cavity function yij(r) = g l] {v)e l} ^ zi{ - r \ 

A similar result holds for the virial route to the equa- 
tion of state, 



Z 



P 



pksT 

= 1 + \p Yl XiX i J dr yv ( r ) r ■ Ve 



-/30. 3 (r) 



, (3.8) 



where P is the pressure. First, note that 
Vfa (r) = 



^4iMe (coBfl ) + %Me(-ooBtf) 



dr dr 
8{z) [fij (r) - tpji (r)]z. 



r 

(3.9) 



Therefore, 
r • Vfa(r) 

and thus 



d ^ (r) e(cos^) + %^e(-cos^) 



dr 



dr 



(3.10) 



ff«(r)^M= I" d(™0) 9ij (r). (3.4) Z=l + \p^ J j dr Vij (r) r^e^ (3AI) 

J-l i,3 
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Finally, let us consider the compressibility route. In a 
mixture, the (dimensionless) isothermal compressibility 
Xt is in general given by 



XT 1 



1 (dp 

k B T \d P/T 



I + h (0) 



(3.12) 



where hij (0) is proportional to the zero wavenumbcr limit 
of the Fourier transform of the total correlation function 
hij(r) = gij(r) - 1, namely 



hij(0) = Py/xixj / dxhij (r) 



= / dv [MO + MO]- (3-13) 

In the specific case of a binary mixture considered here, 
Eq. (|3. 12|) becomes 



_ x 1 + x 2 h n (0) + x 1 h 22 (0) - 2^/xTx 2 'h 12 (0) 
X T = 



i + Mo)J [i + Mo)J -/i? 2 (o) 

(3.14) 

Equations (|3/F]), (j3~TTj) . (|3~T2]) , and (|3~T3)l confirm that 
the knowledge of the two average quantities <7y(r) and 
gji(r) for each pair ij suffices to determine the thermo- 
dynamic quantities. In fact, Eqs. (|3.7[) . p. lip . (|3.12l) . 
and (|3.13p are formally indistinguishable from those cor- 
responding to mixtures with standard isotropic interac- 
tions, except that in our case one generally has (fij(r) 7^ 
ipji(r) and, consequently, g l3 {r) ^g^{r). 

For future convenience, it is expedient to introduce the 
Laplace transform of fg - (r) : 



dr < 



— sr^— 



Its small- s behavior is of the form 2 ^ 

s 2 G lJ ( s ) = l + < S 2 + ffgV + 



where 



H 



(n) 



dr (—r) n rhij(r). 



Thus, Eq. ([3~T3"jl becomes 

ftij(0) = -2irp^/x~iX~j 
Note that Eq. (j3~16)) implies 



H 



(i) 



(i) 



lim s 2 Gij(s) = 1, 

s— >0 



lim 



s 2 %(s)-l 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



IV. THE STICKY LIMIT 

The mapping of the Kern-Frenkel potential with fixed 
patch orientation along the ±z axis onto a binary mixture 
represents a considerable simplification. On the other 
hand, no approximation is involved in this mapping. 

The presence of the original SW interactions for the 
radial part [sec Eq. (|2.3I) ] makes the analytical treatment 
of the problem a formidable task. Progresses can however 
be made by considering the Baxter SHS limit, for which 
a well defined approximate scheme of solution is avail- 
able in the isotropic case for both one-component^ and 
multi-component — ~— fluids. The discussion reported 
below closely follows the analogue for Baxter symmetric 
mixtures! 19 ' 20 

Let us start by rewriting Eq. (|2.6p as 



(4.1) 



00, r < cr, 
Vij{ r ) = { ~ e V' o- <r < Xcr, 
0, r > Act, 



where en = e 22 = £21 = and £12 = e > 0. However, in 
this section we will assume generic energy scales ey. In 
that case, the virial equation of state (|3.1ip becomes 



Z = 1 + 4?7j/((t) - 12//^>"\r,.r,/, 



A 3 -l 



where 77 = 5 pa 3 is the packing fraction, 



(4.2) 
(4.3) 



is the orientational average global cavity function, and 
Uj ee JL ee i (A 3 - 1) (e^ - 1) (4.4) 

is a parameter measuring the degree of "stickiness" of 
the SW interaction <Pij(r). This parameter will be used 
later on to connect results from numerical simulations 
of the actual Janus fluid with analytical results derived 
for asymmetric SHS mixtures. Although Baxter's tem- 
perature parameters r^- are commonly used in the litera- 
ture, we will employ the inverse temperature parameters 
tij = XjYlTij in most of the mathematical expressions. 

In the case of the interaction potential (|4.ip , the energy 
equation of state (|3.7p reduces to 

tiex = -12— ^XjXjejje^ / drr 2 y l} (r). (4.5) 

® - J <7 

h3 



The compressibility equation of state (|3.12p does not sim- 
plify for the SW interaction . 

Since the (orientational average) cavity function (r) 
must be continuous, it follows that 

Vij(r) = Va(r) [e^e(r -a)- (e^« - l) 6(r - Act)] . 

(4.6) 
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Following Baxter's prescription^ we now consider the 
SHS limit 

A -> 1, tij — > oo, tij = — > (A — l)e/ 3e,j = finite, 

(4.7) 

so that the well (|4.1j) becomes infinitely deep and nar- 
row and can be described by a single (inverse) sticki- 
ness parameter nj . Note that in the present Janus case 
(f-ii = £22 = £21 = 0, ex2 = e > 0) one actually has 
tn = h2 — hi = and t\%=t= l/12r. 

In the SHS limit (|I77)l . Eqs. (14~51 . and gU) be- 

come, respectively, 

Z=l + Arjy{a) - 4?y ^ a^iy [3y y (cr) + cry^ (ct)] , 

(4.8) 

u cx = 1 li '/ ^] '',./•; f , , / , , J/, ; • t ; . (4.9) 



Inserting Eq. (|4.14[) into Eq. f|4. 15[) , we obtain the follow- 
ing large- s behavior of Gij(s): 

e as G lJ (s) = ahijy^a) +ay ij (a)s^ 1 

+ [y ij ^) + ay' ij (a)} s~ 2 + G( S " 3 ). 

(4.16) 

A consequence of this is 

lim^oo e as G lj (s) 

lim^K, s [e as Gij (s) - lim s _ i . 00 e a "Gij (s)] y ' 

(4.17) 



V. A HEURISTIC, IMON-PERTURBATIVE ANALYTICAL 
THEORY 

A. A simple approximate scheme within the 
Percus-Yevick closure 



Vijir) = Vij(r) [9(r - a) + ^^+(r - a)] . (4.10) 



In Eq. 



y'iAa) must be interpreted 



which in principle differs from 
limA->i Vij( r )\ _ ^ However, both quantities coin- 
cide in the one-dimensional case2& and are expected to 
coincide in the three-dimensional case as well. This is 
just a consequence of the expected continuity of -£:Vij(r) 
at r = Acr in the SW case.— 

Thermodynamic consistency between the virial and en- 
ergy routes implies 



dur. 



_dZ dZ 



dp d/3 



(4.11) 



Using Eqs. (|4.8p and (|4.9[) and equating the coefficients of 
tij in both sides, the consistency condition f|4. yields 



dr] 



E 

k.C 



X k X£ 



~ tki^- [3y kt (<r) + ayM)} .(4.12) 



For distances r > a, the orientational averages of the 
cavity functions can be Taylor expanded as 

e(r-a)y ij (r) = e(r-a) [y^a) +y' ij (a)(r - a) + ■ ■ ■] . 

(4.13) 

Hence, if we denote by Yy (s) the Laplace transform of 
0(r — er)ry i:/ (r). Eq. (|4.13j) yields for large s 

e^sYijis) = ay l} {a) + [y tJ (a) + ay'^a)] s' 1 + ■ ■ ■ . 

(4.14) 

According to Eqs. (|4.10[) and (|3.15l) . the relationship 
between the Laplace function Gij(s) and Yij(s) is 



Gij(s) = Y i:j {s) + a 1 t tj y l Aa)e 



(4.15) 



The Ornstein-Zernikc (OZ) equation for an anisotropic 
mixture rcadsi^ 

h ij( r ij) = c ij( r ij) + P^x k dr k h ik (v ik )c kj (v kj ) 

k J 

= Cy(ry) + p^x k dr k c ik (r ik )h kj {r kj ), 

k 

(5.1) 

where cy (r^ ) is the direct correlation function. The PY 
closure reads 



1 



(5.2) 



Introducing the averages cf-{r) = Cjj(r) and c^(r) = 
Cji(r) for Cjj(r) in a way similar to Eqs. (|3.3j) and (|3.4[) . 
Eq. (HHJ) yields 



1 



(5.3) 



Thus, the PY closure for the full correlation functions 
dj (r) and gtj (r) translates into an equivalent relation for 
the orientational average functions ctj{r) and <j i3 (r). A 
similar reasoning, on the other hand, is not valid for the 
OZ relation. Multiplying both sides of the first equality 
in Eq. (|5.1[) by (cos fl^) and integrating over cos^- one 
gets 



Mry) = Cii(Hi) + p^x k J dr k 



dr k I d (cos Qij) 
o 



xh ik (rik)ckj(r k j). 



(5.4) 



The same result is obtained if we start from the sec- 
ond equality in Eq. (|5.ip . multiply by 0(— cos %), in- 
tegrate over cos 8ij, and make the changes — > — r 



r ik —> —r ik , and i j. Equation (|5 .4[) shows that in the 
case of anisotropic potentials of the form (|2.5[) the OZ 
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equation does not reduce to a closed equation involving 
the averages hij(rij) and CijiTij) only, as remarked. 

In order to obtain a closed theory, we adopt the heuris- 
tic mean-field decoupling approximation 



d (cos fly) h ik (r ik )c kj {r kj ) -> 



dr k h lk (r ik ) c kj (r kj ) . 



(5.5) 



Under these conditions, the true OZ relation (|5.4p is re- 
placed by the pseudo-OZ relation 



hjjjrjj) = Cjjjrij) + P/ j dY k h lk (r ik )c kj (r k] ). 

k ■* 

(5.6) 

This can then be closed by the PY equation (|5 .3[) and 
standard theory applies. An alternative and equivalent 
view is to consider Cy(r) not as the orientational aver- 
age of the true direct correlation function Cy (r) but as 
exactly defined by Eq. (|5.6[) . Within this interpretation, 
Eq. (|5.3|) then represents a pscudo-PY closure not deriv- 
able from the true PY closure (|5.2I) . 

Within the above interpretation, it is then important 
to bear in mind that the functions g^ (r) obtained from 
the solution of a combination of Eqs. (|5.3[) and (|5.6[) 
are not equivalent to the orientational averages of the 
functions <?jj(r) obtained from the solution of the true 
PY problem posed by Eqs. (|5.ip and complemented by 
the PY condition (|5.2|) . As a consequence, the solutions 
to Eqs. (|5.3[) and (|5.6|) are not expected to provide the 
exact g t j (r) to first order in p, in contrast to the true 
PY problem. This is an interesting nuance that will be 
further discussed in the next sections. 

The main advantage of the approximate OZ relation 
(|5.6[) in the case of anisotropic interactions of the form 
(|2.5[) is that it allows to transform the obtention of an 
anisotropic function gij{v), but symmetric in the sense 
gij(r) = r), into the obtention of an isotropic func- 
tion g~ij{r), but asymmetric since g~ij(r) ^ 9ji(r). In the 
case of the anisotropic SHS potential defined above, we 
can exploit the known solution of the PY equation for 
isotropic SHS mixtures to construct the solution of the 
set made of Eqs. (|5.3[) and (|5.6[) . This is done in the next 
subsection by following the RFA methodology. 



B. RFA method for SHS 

Henceforth, for the sake of simplicity, we take a = 1 
as length unit. The aim of this section is to extend 
the RFA approximation proposed for symmetric SHS 
mixture o ' 25 to the asymmetric case. 

We start with the one-component case^ Let us intro- 
duce an auxiliary function F(s) related to the Laplace 
transform G(s) of rg(r) by 



G(s) 



1 



2tt F{s) + pe- 



(5.7) 



Next step is to approximate F(s) by a rational function, 

S(s) 



F(s) 



L( S y 



with S(s) = + SWfl + SW s 2 + s 3 and 
L(s)=lW+lW S + lWs 2 . 



(5.8) 



(5.9) 



Note that lim^oc F(s)/s = = finite, so that 

linis-yco e s G(s) = finite, in agreement with Eq. (|4.1 61) . 
Furthermore, Eq. (|3.19[) requires F(s) + pe~ s = 0(s 3 ), 
so that S<°) = - P L(°), S« = p(L^-L^), = 
p (iW - - L^). Taking all of this into account, 

Eq. (|5.7p can be rewritten as 



G(s) = 



L(s) 



2ns 2 1 - A(s) ' 



(5.10) 



where 



Ms) = 4 



(l-e- s )L(s)-L^s+ Ql^-LW 



(5.11) 

In the case of a mixture, G(s), L(s), and A(s) become 
matrices and Eq. (|5.10[) is generalized as 

G^( S ) = ^(l( S ).[I-A( S )]- 1 ).. j (5.12) 
where I is the identity matrix and 



Lv{8) = LM+L$>a + L$>*, 



(5.13) 



A ij (s)=p^ (1 



Ir(°) t 



L i:j (s) - L 
(i 



(5.14) 



Note that lim^o Aij (s) = finite, so that 
lim s _^o s 2 Gij{s) = finite 7^ by construction. Analo- 
gously, lim s _ >00 e s Gij(s) = finite also by construction. 

The coefficients Lyy, L^]\ and are determined by 
i j i j *j ° 

enforcing the exact conditions (|3.19[) . (|3.20[) , and (|4.17l) . 

The details of the derivation are presented in Appendix 
|B]and here we present the final results. The coefficients 
and do not depend upon the first index i and 
can be expressed as linear functions of the coefficients 

{L k % 



(o) 



2tt 



1 + 2?; 12?7 



(1-1,)= 



L 



8 



and the coefficients L$ obey the ciosed set of quadratic Next, Eqs. (|5.15j) and (|5.16|) yield 



equations 



L 



12) 



= 2tt 



1 + r]/2 677 



n 

(2) r (2) 



* fc (^ + 4?) 



(5.17) 



This closes the problem. Once L\- arc known, the con- 
tact values are given by 



Ml) 



I 



(2) 



(5.18) 



Although here we have taken into account that all the 
diameters are equal (er^ = a = 1), the above scheme can 
be easily generalized to the case of different diameters 
with the additive rule a\j = (04 + o~j)/2. For symmetric 
interactions (i.e., Uj = tji) one recovers the PY solution 
of SHS mixtures for any number of components. 22 ^ 4 It is 
shown in Appendix [C] that the pair correlation functions 

(r) derived here are indeed the solution to the PY-like 
problem posed by Eqs. (|5.3[) and (|5.6[) . 



C. Case of interest: tu = £22 = £21 = 



In the general scheme described by Eqs. (|5. 1 2[) — (|5. 18[) . 
four different stickiness parameters (tu, ti2, t2i, and £22) 
are in principle possible. With the convention that in 
tij particle i is always located below particle j, we might 
consider the simplest possibility of having only one SHS 
interaction ti2 = t = l/12r and all other HS interactions 
(til = *22 = t-2i = 0), as illustrated in Fig. [2j This is 
clearly an intermediate case between a full SHS model 
(t^ = t = l/12r) and a full HS model (tij = 0), with a 
predominance of repulsive HS interactions with respect 
to attractive SHS interactions. This is meant to model 
the intermediate nature of the original anisotropic Kern- 
Frenkel potential that interpolates between a SW and 
a HS isotropic potentials upon decreasing the coverage, 
that is, the fraction of the SW surface patch with respect 
to the full surface of the sphere. 



r(0) 

2tt 



L 



(0) 



21 
2tt 



1 +27? 

(1 - 7?) 2 



(5.21) 



r(0) r(0) 



J 12 

2tt 



1 + 2r7 12r)t 



2tt 



1 



17 



Ziy 12 (l), (5.22) 



r(l) 

2tt 



(i) 



21 
2ir 



(I-77) 2 ' 



J \2 

2tt 



22 

2tt 



1 + g/g 

(I-77) 2 1-7? 



W - ,,, 
£12/12(1)- 



(5.23) 



(5.24) 



Once the functions Lij(s) are fully determined, Eq. 
(|5.12j) provides the Laplace transforms G^ (s). From Eq. 
(I4T51) it follows that F u (s) = G*n(s), r 22 (s) = G 22 (s), 
*2i(s) = G 2 i(s), and 



r 12 («) = Gia(«)-*y 12 (l)e- 



(5.25) 



A numerical inverse Laplace transform 2 ^ allows one to 

obtain <? u (r), 7722M, ff2i0): and yi2( r ) for an Y pack- 
ing fraction 7?, stickiness parameter t = l/12r, and mole 
fractional. In what follows, we will omit explicit expres- 
sions related to g 2 2(?") since they can be readily obtained 
from 3n(r) by the exchange x\ O x 2 . 

The contact values g i? -(l + ) — |fa (l) with (£, j) 7^ (1, 2) 
cannot be obtained from Eq. (|5.18j) , unless the associated 
tij are first assumed to be nonzero and then the limit 
tij — > is taken. A more direct method is to realize that, 
if = 0, Eq. (|4TTd| gives 

&,(!+) = lim e s sGij(s), ^ (1,2). (5.26) 



The results are 



5ii(i + ) = yn(i) = 



(1 - nf 



52^-1/12(1), (5.27) 



ff 2 i(i + ) = y 2 i(i) = 



T — 77 1 



1 + ^/2 
(I-7?) 2 



(5.28) 



1. Structural properties 



If *ii = £22 = hi = 0, Eq. (I5~T71) implies L$ = L 22 > = 



(2) 



J 2i =0. As a consequence, Eq. (|5.17p for i 

(2) 



1 and 



j = 2 yields a linear equation for L 12 whose solution is 

(5.19) 



1 — 77 1 — 7? + 6r/t 
According to Eq. (|5.18|) , 

y (1) = 1 + 71/2 (l - 67lt 

(1 — 7/) 2 \ 1 — 7/ + 6rjt 



(5.20) 



17(1) = 



l + *7/2 
(I-77) 2 



1 



2:12:2 



12??t 



1 — 77 + 67?t 



(5.29) 



It is interesting to note the property 5n(l + ) + g 2 2(l + ) = 

yi 2 (l)+5 2 i(l + )- 

To obtain the equation of state from the virial route 
we will need the derivative j7i 2 (l)- Expanding e s G\2(s) 
in powers of s _1 and using Eq. ()4.16|) . one gets 



gjgjl) 
¥12(1) 



(1 - 7?) 2 
9 1- 



3t 



2 - 4?; - 7r/ 2 
1+77/2 



12a;iX 2 77 



2 1 



2 1 + ?7/2 



(5.30) 



2. Thermodynamic properties 



Virial route. According to Eq. (j4,8|) . 
Z v = l + 4r ] y(l)-4x 1 X2Vt [3y 12 (l) + y' 12 (l)] 



^HS 



1 



%i 2 (i) + y'i 2 (i) 



(5.31) 

where the superscript v denotes the virial route and 

1 + 2rj + 3r ? 2 



(5.32) 



that, for an equimolar mixture (x\ — X2 = h), reduces 



to 



Xt 



[(1 - 77) 2 (1 + 2/ ? ) + 3 ?? i (2 + 5r ? - fr? 2 )]' 



(1 - 77)5(1 - r/ + 6nt) [(1 - 7 ? ) 2 + 3^(4 - t))] ' 

(5.38) 

The associated compressibility factor is then 



Z c = ~ P 'drfxT^rf). 
V Jo 



(5.39) 



The above integral has an analytical solution, but it is 
too cumbersome to be displayed here. 



is the compressibility factor predicted by the virial route 
in the PY approximation. 

b. Energy route. From Eq. (|4.9p we have 



— = -12x 1 x 2 r]ty 12 (l). 
e 



(5.33) 



The compressibility factor can be obtained from u ox via 
the thermodynamic relation (|4. 1 1 [) . which in our case 
reads 



d u cx /e 
dr] 



1 dZ _ dZ 
eif ~ t ~dt' 



(5.34) 



Thus, the compressibility factor derived from the energy 
route is 



Z" 



^HS 



7 



HS 



d_ 

dr] 



3xix 2 - 



* dt , u ex (y,t')/e 



In 1 



4^i 2 (l) + 



6yt 
1-') 



(5.35) 



where Z^ s plays the role of an integration constant and 
thus it can be chosen arbitrarily. It can be show n 29 i 30 
that the energy and the virial routes coincide when the 
HS system is the limit of a square-shoulder interac- 
tion with vanishing shoulder width. From that point of 
view one should take Z^ s = Z^ s in Eq. (|5.35[) . On the 
other hand, a better description is expected from the 
Carnahan-Starling (CS) equation of state 



7CS _ 

^HS — 



1 



+ 1] + 1]" 



(1-77)3 



(5.36) 



3. Low-density expansion 

In the standard case of SHS mixtures with symmetric 
coefficients in the potential parameters, the PY closure 
is known to reproduce the exact cavity functions to first 
order in density and thus the third virial coefficient (sec 
Appendix I A 21) . However, this needs not be the case in 
the present asymmetric case, as further discussed be- 
low. Note that here, "exact" still refers to the simplified 
problem (orientational average+sticky limit) of Sections 
|mi and EI 

The expansion to first order in density of the Laplace 
transforms Yij(s) obtained from Eqs. (|4.15[) . (|5,12|> 
(|5H|) . and ([5T9"1) - (|5~^4]) is 



(5.40) 



where the expressions of the first-order coefficients 
Yu' (s) will be omitted here. Laplace inversion yields 



_(l), \ X 



(5.41) 



where y^\r) 



are the exact first-order functions 



given by Eqs. (|A37j> — (|A39|) and the deviations Ay£\ 



arc 



Ay$(r)=Q(2-r)x 



2f 



2 — cos -, 
r 2 



(5.42) 



Ay$(r) = B(2-r)t fey/l - r 2 /4 - rcos" 1 , (5.43) 



Henceforth we will take Z^ s 



c. 



Compressibility route. Expanding s Gy(s) in 



Ay 2 \\r) = -Ay£(r) 



(5.44) 



powers of s it is straightforward to obtain h[] ] from Eq. In thc case of the g lobal quantity y W(r) the result is 



([3"36]) . This allows one to use Eqs. (pTH)) and (f3TT5|) to 

get the inverse susceptibility as 



y [1 \r) = t L, (r) 



(5.45) 



X T 



l + 2i] l + 2r)-24x 1 x 2 tr](l-ri)y 12 (l) 



(1 - rj)< 



1 - x±x 2 



12t V (l+ n /2)y 12 (l) 



1+277+36x1X2 tr) 2 y 12 (l) 



(5.37) 



where ^ 1 ^ C 7 ^) I i s given by Eq. (|A"4"0|) and 

Ay^'M = Q(2-r)x 1 x 2 — cos" 1 -. 

r 2 



(5.46) 
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While the main qualitative features of the exact cav- 
ity function are preserved, there exist quantitative dif- 
ferences. The first-order functions Vyiir), 2/22* ( r )> anc ^ 
y( x )(r) predicted by the RFA account for the exact co- 
efficient of t but do not include the exact term of or- 
der t 2 proportional to r _1 cos -1 (r/2). In the case of 
y^2 ( r ) and y^ (r) the exact term of order t propor- 
tional to 2-y/l — r 2 /4 — rcos _1 (r/2) is lacking. Also, 

while the combination y^ (r) + ( r ) — V12 ( r ) ~ V21 ( r ) 
vanishes in the RFA, the exact result is proportional to 
£2 r -i C os _1 (r/2). In short, RFA correctly accounts for 



TABLE I. Location of the critical point in the RFA, according 
to different routes. 



Route 




Vc 


z c 


virial, Eq. 1)5.3] \ 


0.02050 


0.1941 


0.3685 


energy, Eq. (|5.35|) 


0.0008606 


0.2779 


0.2906 


hybrid virial-enerev. Ea. (15.501) 


0.01504 


0.1878 


0.3441 



the polynomial terms in yiV (r) 



Within the PY approximation, the critical behavior 
of the original one-component Baxter SHS model was 

15 



but misses the non- studied using the compressibility and virial routes,— as 



polynomial terms. 

As for the thermodynamic quantities, expansion of 
Eqs. ([OTj) . ([535] . and (jO^ gives 



Z" = 1 + 4(1 - ix 1 x 2 t)r s 



10 



1 — 6x\X2t 1 — 



4 



(5.47) 



Z u = 1 + 4(1 - 'ix l x 2 t)r) + 10 

+o{rf), 

Z c = 1 + 4(1 - 3x!X 2 t) 7] + 10 
+0{rf). 



6 

1 — OX\X2t 1 — — t 



(5.48) 



1 — %x\x 2 t [ 1 — —t 



n 

(5.49) 



Comparison with the exact third virial coefficient, Eq. 
(|A50|) . shows that the coefficient of t 2 is not correct, 
with the exact factor 4 — 3%/3/tt — 2.35 replaced by 2, 
3, and 4 in Eqs. (|5.47[) - (|5.49[) . respectively. One con- 
sequence is that the virial and energy routes predict the 
third virial coefficient much better than the compress- 
ibility route. A possible improvement is through the 
interpolation formula 



Z v - U = a (Z v + Zg| - Z^ s ) + (1 - a)Z u , (5.50) 



where a = 3\/3/ir — 1 ~ 0.65 with the proviso that 
Z£ s = Zgf in Eq. (j5.35|) . Equation (|5.50j) then reduces 
to the CS equation of state if t = and reproduces the 
exact third virial coefficient when t 7^ 0. 



4. Phase transition and critical point 

In the limit of isotropic interaction (t^j = t), our 
model reduces to the usual SHS Baxter adhesive one- 
component model. In spite of the fact that the model 
is, strictly speaking, known to be pathological^ it dis- 
plays a critical behavior that was numerically studied in 
some details by MC techniques . 32 ' 33 The corresponding 
binary mixture also displays well defined critical proper- 
ties that, interestingly, are even immune by any patho- 
logical behavior^ 



well as the energy routed in the latter case with the 
implicit assumption Z£ s = Z§§. Numerical simulations 
indicate that the critical point found through the energy 
route is the closest to numerical simulation results i 32 ' 33 
As the present specific model ( with, t%j = t5n8j2) is, in 
some sense, intermediate between the fully isotropic Bax- 
ter SHS one-component model (that has a full, albeit pe- 
culiar, gas-liquid transition) and the equally isotropic HS 
model (that, lacking any attractive part in the potential, 
cannot have any gas- liquid transition) , it is then interest- 
ing to ask whether in the equimolar case {x\ — x 2 = ^) 
it still presents a critical gas-liquid transition. 

The answer depends on the route followed to obtain the 
pressure. As seen from Eq. (|5.38[) , the compressibility 
route yields a positive definite x^ 1 , so that no critical 
point is predicted by this route. On the other hand, 
an analysis of the virial [Eq. (|5.31[) ]. energy [Eq. (|5.35[) 
with Z£ s = and hybrid virial-energy [Eq. (|5.50[) ] 

equations of state reveals the existence of van der Waals 
loops with the respective critical points shown in Table 
HI The energy route predicts a critical value r c about half 
the values predicted by the other two routes. 

As an illustration, Fig. 2] shows the binodal and a few 
isotherms, as obtained from the virial route. 



5. A modified approximation 

The failure of the RFA to reproduce the exact cavity 
functions to first order in density (and hence the third 
virial coefficient) for asymmetric interactions (t^ =/= tji) 
reveals the price paid for using the orientationally aver- 
aged quantities (r) instead of the true pair correlation 
functions <?y (r). 

A simple way of getting around this drawback for suf- 
ficiently low values of both r\ and t consists of modifying 
the RFA as follows: 



(5.51) 



where the functions Aj/^^r) are given by Eqs. (|5.42[) - 
(|5.44[) . We will refer to this as the modified rational- 
function approximation (mRFA). Note that Eq. (|5.51[) 

implies that g n (r) -> g i:j (r) + Hyf) {r)p, except if (i, j) = 
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FIG. 4. Binodals from the RFA virial route in the equimolar 
xi = \ case. The phase diagram is depicted in the (77, r) plane 
(top panel) and in the (r\~ x ,r]Z v ) plane (bottom panel). The 
critical point is found at 77c ~ 0.1941, r c ~ 0.02050. 



(1,2), in which case g 12 (r) -4 g 12 (r) + Ay\ 2 '(r)p + 
Ay^(l)S+(r-l)pt. 

Since the extra terms in Eq. (|5.51[) are proportional to 
t or t 2 , this modification can produce poor results for 
sufficiently large stickiness (say, t > 1) as, for instance, 
near the critical point. 



VI. NUMERICAL CALCULATIONS 

A. Details of the simulations 

In order to check the theoretical predictions previously 
reported, we have performed NVT (isochoric-isothermal) 
MC simulations using the Kern-Frenkel potential defined 
in Eqs. (|2. ![) — (|2.4[) with a single attractive SW patch 
(green in the color code of Fig. [T]) covering one of the 
two hemispheres, and with up-down symmetry as de- 
picted in Fig. [5J Particles are then not allowed to rotate 
around but only to translate rigidly. 

The model is completely defined by specifying the rel- 



ative width A — 1, the concentration of one species (mole 
fraction) x\ = 1 — X2, the reduced density p* = pa 3 , and 
the reduced temperature T* = kBT/e. 

In order to make sensible comparison with the RFA 
theoretical predictions, we have selected the value A — 1 = 
0.05 as a well width, which is known to be well repre- 
sented by the SHS limit^i and use Baxter's temperature 

parameter r = [4(A 3 - 1) (e 1/T * - l)] ~* [sec Eq. (|Q|) ] 
instead of T* . 

During the simulations we have computed the oricn- 
tational averaged pair correlation functions defined by 
Eqs. (|3.3[) and (|3.4[) . accumulating separate histograms 
when Z2 — z\ > or z\ — z 2 > in order to distinguish 
between functions = 5i2( r ) an d 52i( r ) = Si^M- 

The compressibility factor Z = f3P/ p has been evalu- 
ated from the values of TJijir) at r = a and r = \a by 
following Eq. (|4.2[) with = (12T)~ 1 5n5j2, and the re- 
duced excess internal energy per particle u* x = u cx /e has 
been evaluated directly from simulations. 

For most simulations, we used N = 500 particles, 
periodic boundary conditions, an equilibration time of 
around 10 5 MC steps, and a production time of about 
10 6 MC steps, with a maximum displacement of particles 
selected as to maintain an acceptance ratio to about 0.5. 

B. Results for non-equimolar binary mixtures 

As a preliminary attempt, we consider a binary mix- 
ture under non-equimolar conditions, to avoid possible 
pathologies arising from the symmetry of the two com- 
ponents akin to those occurring in ionic systems. As we 
shall see below, no such pathologies are found. 

In the present case, we consider a system with xi = 1/5 
and X2 = 1—xi = 4/5, so that the majority of the spheres 
have (green) attractive patches pointing (— z). 

A snapshot of an equilibrated configuration is shown in 
Fig. O This configuration was obtained using N = 500 
particles at p* = 0.3 and Baxter temperature r = 0.1 
(corresponding to T* ~ 0.354). Note that the above 
chosen state point (p* = 0.3 and r = 0.1) lies well inside 
the critical region of the full Baxter SHS adhesive model 
as obtained from direct MC simulations ; 32 ' 33 although of 
course the present case is expected to display a different 
behavior as only a fraction of about X1X2 = 4/25 of the 
pair contacts are attractive. 

A good insight on the structural properties of the 
system can be obtained from the computation of the 
radial distribution functions j u (r). ffi2( r ) ~ 5i2( r )i 
512M = 521 ( r )> and 522 M- This is reported in Fig. [6] 
for a state point at density p* = 0.5 and Baxter temper- 
ature t = 0.2 (corresponding to T* ~ 0.457). Note that 
in the case of the pair (1,2) what is actually plotted is the 
cavity function y 12 (r) rather than 3 12 (^), as explained in 
the caption of Fig. [51 

The relatively low value r = 0.2 gives rise to clearly 
distinct features of the four MC functions ~g~ij(r) (which 
would collapse to a common HS distribution function 
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FIG. 5. Snapshot of an equilibrated MC simulation under 
non-equimolar conditions (xi = 1/5) with Baxter tempera- 
ture t = 0.1 and density p* = 0.3. In the simulations we used 
a total number of particles N = 500. 



in the high-tempcraturc limit r — > oo). Wc observe 
that g 22 (r) ~ g 21 (r) > > V12 ( r ) in the region 

1 < r ^ 1-5. Moreover, g n (r) and <7i2(?") exhibit a rapid 
change around r = 2. This is because when a pair (1,1) 
is separated a distance r ~ 2 there is enough room to fit 
a particle of species 2 in between and that particle will 
interact attractively with the particle of the pair (1,1) be- 
low it. In the case of the pair (1,2) separated a distance 
r 2, the intermediate particle can be cither of species 
1 (interacting attractively with the particle of species 2 
above it) or of species 2 (interacting attractively with the 
particle of species 1 below it). The same argument ap- 
plies to a pair (2, 2) separated a distance r ~ 2, but in 
that case the intermediate particle must be of species 1 
to produce an attractive interaction; since the concentra- 
tion of species 1 is four times smaller than that of species 
2, the rapid change of g 22 (r) around r = 2 is much less 
apparent than that of j u (r) and 312(f) in Fig. [6] On the 
other hand, in a pair (2, 1) separated a distance r w 2 
an intermediate particle of either species 1 or of species 2 
does not create any attraction and thus <? 2 i ( r ) i s rather 
smooth at r = 2. In short, the pair correlation function 
g 2 i(r) exhibits HS-likc features, g 12 (r) exhibits SW-likc 
features (very high values in the region 1 < r < A and 
discontinuity at r = A due to the direct SW interaction; 
rapid change around r = 2 due to indirect SW interac- 
tion), while ~g~i\(r) and g 22 (r) exhibit intermediate fea- 
tures (rapid change around r = 2 due to indirect SW 
interaction). 

It is rewarding to notice how well the MC results 
are reproduced at a semi-quantitative level by the RFA 
theory (top panel of Fig. in spite of the various ap- 
proximations involved. In this respect, it is worth recall- 
ing that while MC simulations deal with the real Kern- 
Frcnkcl potential, albeit with constrained angular orien- 
tations, the RFA theory deals with the asymmetric bi- 
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FIG. 6. Comparison between MC simulations and the the- 
oretical predictions from RFA (top) and mRFA (bottom) 
for the orientational averaged distribution functions j u (r), 
y 12 (r), <7 2 i( r )> arm 922 ( r ) under non-equimolar conditions 
(xi = 1/5) at density p* = 0.5 and Baxter temperature 
r = 0.2. The dashed vertical line indicates the range r = 
A = 1.05 of the (1, 2) SW interaction used in the simula- 
tions. Note that the radial distribution function (? 12 (r) is 
obtained in the MC case by multiplying y 12 (r) in the region 
1 < r < A = 1.05 by the factor e 1/T * ~ 8.93; in the theoreti- 
cal cases (SHS limit) g 12 (r) is obtained by adding the singular 
term {12r)- l y X2 {l)8+{r - 1) to y 12 (r). 



nary mixture resulting from the mapping described in 
Section mi and this represents an indirect test of the cor- 
rectness of the procedure. In addition, the RFA does 
not attempt to describe the true SW interaction (i.e., fi- 
nite A - 1 and T*) but the SHS limit (A - 1 -> and 
T* — > with finite r). This limit replaces the high jump 
of #12 (f) in the region 1 < r < A by a Dirac's delta at 
r = 1 + and the rapid change of g 12(f), gn(r), and g 2 2( r ) 
around r = 2 by a kink. Finally, the RFA worked out in 
Sec. IV Bl results from a heuristic generalization to asym- 
metric mixtures (r^ 7^ Tji) of the PY exact solution for 
SHS symmetric mixtures (t^ = t^)^" -22 ' 24 but it is not 
the solution of the PY theory for the asymmetric prob- 
lem, as discussed in Sec. IV Al As a matter of fact, the top 
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panel of Fig. [6] shows that some of the drawbacks of the 
RFA observed to first order in density in Sec. IV C 31 [see 
Eqs. (|5.41j) - (|5.44j) ] remain at finite density: in the region 
1 < r S 1-5 the RFA underestimates 1/12(1*), ~g~ii{ r )i an< i 
g 22 (?"), while it overestimates g 21 (r). These discrepancies 
are widely remedied, at least in the region 1 < r < 1.25, 
by the mRFA approach [see Eq. (|5.51[) ]. as shown in the 
bottom panel of Fig. [SJ In particular, the contact val- 
ues are well accounted for by the mRFA, as well as the 
property <? 22 (?~) — ff2i( r )- We have observed that the 
limitations of the correlation functions Tj^ (r) predicted 
by the RFA become more important as the density and, 
especially, the stickiness increase and in those cases the 
mRFA version does not help much since the correction 
terms Ay[^\r)p, being proportional to p and to t or t 2 , 
become too large. 
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FIG. 7. Comparison of MC simulations and RFA theory for 
the thermodynamics. Both the compressibility factor Z — 
/3P/p (top) and the excess internal energy per particle it ox /e 
(bottom) are displayed as functions of density for the non- 
equimolar case xi = 1/5 and for Baxter temperature r = 0.1. 
In the case of the compressibility factor (top), results for all 
four routes (compressibility, virial, energy, and hybrid virial- 
energy) are reported. 

Next we consider thermodynamic quantities, as repre- 
sented by the compressibility factor Z = j3P/ p and the 



excess internal energy per particle u cx /e, both directly 
accessible from NVT numerical MC simulations. These 
quantities are depicted in Fig. [7J as functions of the re- 
duced density p* and for a Baxter temperature r = 0.1. 
In both cases, the results for the RFA theory are also in- 
cluded. In the case of the compressibility factor, all four 
routes are displayed: compressibility [Eqs. (|5.20[) , (|5.37l) . 
and flOg) )], vir ial [ Eqs. (E TUI), (|QTj)) . and (jOTj) ]. en- 
ergy [Eq. (ph2"0)) and ([Q5)l with Z£ s = Zgf], and hybrid 
virial-energy [Eq. (|5.50|) ]. In the case of u ex /e, only the 
genuine energy route, Eq. (|5.33|) . is considered. Note 
that all RFA thermodynamic quantities, including Eq. 
(|5.39[) . have explicit analytical expressions. 

The top panel of Fig. [7| shows that up to p* « 0.7 the 
MC data for the compressibility factor are well predicted 
by the theoretical Z v and, especially Z u and Z v - U . Be- 
yond that point, the numerical results are bracketed by 
the compressibility route, that overestimates the pres- 
sure, and the hybrid virial-energy route, that on the 
contrary underestimates it. It is interesting to note 
that, while Z v < Z v ' u < Z u to second order in den- 
sity [cf. Eqs. (pT4"?]) . (pT^ . and (j5"30"]) ]. the difference 
Z v — grows with density more rapidly than the dif- 
ference Z u — Z^ s and so both quantities cross at a certain 
density (p* ~ 0.567 if x x = 1/5 and r = 0.1). There- 
fore, even though Z v < Z u (because Z£ s < Zgf), Z v ' u 
is no longer bracketed by Z v and Z u beyond that den- 
sity (p* ~ 0.567 in the case of Fig. [7]). On balance, the 
virial-energy route appears to be the most effective one 
in reproducing the numerical simulations results of the 
pressure at x\ = 1/5 and r = 0.1. 

As for the internal energy, the bottom panel of Fig. [7] 
shows that the RFA underestimates its magnitude as a 
direct consequence of the underestimate of the contact 
value y 12 (l) [see Eq. (|5.33|) ]. Although not shown in Fig. 
[Jj we have checked that the internal energy per particle 
obtained from the virial equation of state (|5.31[) via the 
thermodynamic relation (|5.34[) exhibits a better agree- 
ment with the simulation data than the direct energy 
route. 



C. Results for equimolar binary mixtures 

Having rationalized the non-equimolar case, the 
equimolar (x\ = x 2 = 1/2) case can now be safely tack- 
led. The equimolarity condition makes the system be 
more akin to the original Janus model (see Fig. [JJ since 
both spin orientations are equally represented. 

We start with the snapshot of an equilibrated con- 
figuration at density p* = 0.3 and Baxter temperature 
t = 0.1, that are the same values used in the non- 
equimolar case. It can be visually inspected that, in 
contrast to the non-equimolar case of Fig. [SJ the num- 
ber of particles with spin up matches that with spin 
down. This equimolar condition then facilitates the in- 
terpretation of the corresponding structural properties, 
as illustrated by the radial distribution function (r) 
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FIG. 8. Same as in Fig. [5] but for an equimolar mixture 
{xi =x 2 = 1/2). 



given in Fig. [9] 

This was obtained at a Baxter temperature r = 0.2 
and a density p* = 0.5, a state point that is expected 
to be outside the coexistence curve (see below), but in- 
side the liquid region. Again, this is the same state 
point as the non-equimolar case previously discussed. 
Now j u (r) = g 22 (r) (independently computed) as it 
should. Notice that the main features commented before 
in connection with Fig. [5] persist. In particular, g 21 (r) > 
Vn( r ) =j22( r ) > Vi2( r ) in the region 1 < r < 1.5, 
5n( r ) = 522 ( r ) ancl 3i2( r ) present rapid changes around 
r = 2, and g 21 ( r ) exhibits a HS-like shape. Also, as be- 
fore, the RFA captures quite well the behaviors of the 
correlation functions (especially noteworthy in the case 
of g xl = 1)22) ■ On the other hand, the RFA tends to 
underestimate y 12 (r) and g~n(r) = 522 ( r ) an d to overes- 
timate 5 21 (r) in the region 1 < r < 1.5. The use of the 
modified version (mREA) partially corrects those discrep- 
ancies near contact, although the general behavior only 
improves in the case of g 21 (r). 

Comparison between Figs. |5] and |H] shows that y^ir) 
and g 21 (r) are very weakly affected by the change in com- 
position. In fact, the spatial correlations between parti- 
cles of species 1 and 2 mediated by a third particle (i.e., 
to first order in density) depends strongly on which par- 
ticle (1 or 2) is above o below the other one but not on 
the nature of the third intermediate particle, as made ex- 
plicit by Eqs. (jA38]) and (|A39|) . Of course, higher-order 
terms (i.e., two or more intermediate particles) create a 
composition-dependence on y 12 (r) and g 2 i{ r )i but this 
effect seems to be rather weak. On the contrary, the mi- 
nority pair increases its correlation function 3 n (r), while 
the majority pair decreases its correlation function g 22 (r) 
in the region 1 < r < 1.5 when the composition becomes 
more balanced. Again, this can be qualitatively under- 
stood by the exact results to first order in density [see 

Eq. 4M3)]. 
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FIG. 9. Same as in Fig. [6l but for an equimolar mixture 
{xi = x 2 = 1/2). 



D. Preliminary results on the critical behavior 

One of the most interesting and intriguing predictions 
of the RFA is the existence of a gas-liquid transition in 
the equimolar model, despite the fact that only one of the 
four classes of interactions is attractive (sec Sec. IVCif . 
The clusivcncss of this prediction is reflected by the fact 
that the compressibility route does not account for a crit- 
ical point and, although the virial and energy routes do, 
they widely differ in its location, as seen in Table HI In 
this region of very low values of r the hybrid virial-energy 
equation of state is dominated by the virial one and thus 
the corresponding critical point is not far from the virial 
value. 

A simple heuristic argument can be used to support 
the existence of a true critical point in our model. Ac- 
cording to the Noro-Frenkel (NF) generalized principle of 
corresponding states,— the critical temperatures of dif- 
ferent systems of particles interacting via a hard-core po- 
tential plus a short-range attraction are such that the 
reduced second virial coefficient i? 2 = £? 2 /£?H s has a 
common value £? 2 C ~ — 1.21. In our model, the reduced 
second virial coefficient is _B 2 = 1 — 3</4 = 1 — l/16r [see 
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FIG. 10. MC simulation data for the scaled pressure r\Z — 
fo- 3 /3P as a function of p* at r = 0.030 (top panel), 0.0205 
(middle panel), and r = 0.018 (bottom panel). Densities 
higher than p* = 0.7 are not shown because at these very low 
temperatures the particles tend to overlap their SW shells and 
then the calculations slow down considerably. Also shown are 
the theoretical results for the four routes of the RFA. 



Eq. (|A49[) ] . Thus, assuming the NF ansatz, the critical 
point would correspond to t^ f ~ 0.028, a value higher 
than but comparable to that listed in Table U from the 
virial route. 

From the computational point of view, a direct as- 
sessment on the existence of a gas-liquid transition in 
the present model is not a straightforward task. Unlike 



the original SHS Baxter model, a Gibbs ensemble MC 
(GEMC) calculation for a binary mixture is required to 
find the coexistence lines. We are currently pursuing this 
analysis that will be reported elsewhere. As a very pre- 
liminary study, we here report NVT results (this time 
with N — 5000 statistics accumulated over 5 x 10 8 MC 
steps) with values of the Baxter temperature close to 
the critical value r^ F ~ 0.028 expected on the basis 
of the NF conjecture. More specifically, we have con- 
sidered t = 0.030, 0.0205, and 0.018 (corresponding to 
T* ~ 0.251, 0.229, and 0.223, respectively). The numeri- 
cal results for the pressure, along with the RFA theoret- 
ical predictions, are displayed in Fig. [TU] 

We observe that at r = 0.030 (top panel) the four 
theoretical routes clearly indicate a single-phase gas-like 
behavior with a monotonic increase of the pressure as 
a function of the density, in consistence with the value 
t c ~ 0.0205 obtained from the RFA virial route. On 
the other hand, the MC data show a practically constant 
pressure between p* = 0.2 and p* = 0.4, which is sug- 
gestive of r = 0.030 being close to the critical isotherm 
(remember that rf F ~ 0.028). The middle panel has 
been chosen to represent the critical isotherm predicted 
by the RFA-virial equation of state. In that case, the 
simulation data present a clear van der Waals loop with 
even negative pressures around the minimum. A simi- 
lar behavior is observed at r = 0.018 (bottom panel), 
except that now the RFA-virial isotherm also presents a 
visible van der Waals loop. Whereas the observation of 
negative values of isothermal compressibility in the MC 
simulations can be attributed to finite-size effects and are 
expected to disappear in the thermodynamic limit, these 
preliminary results are highly supportive of the existence 
of a gas-liquid phase transition in our model with a crit- 
ical (Baxter) temperature r c « 0.03. 



VII. CONCLUSIONS AND FUTURE PERSPECTIVES 

We have studied thermophysical and structural prop- 
erties of a Janus fluid having constrained orientations 
for the attractive hemisphere. The Janus fluid has been 
modeled using a Kern-Frenkel potential with a single SW 
patch pointing either up or down, and studied using nu- 
merical NVT MC simulations. 

The above model has been mapped onto an asymmet- 
ric binary mixture where the only memory of the origi- 
nal anisotropic potential stems from the relative position 
along the z-axis of particles of the two species 1 and 2. In 
this way, only one [(1,2) with our choice of labels] out of 
the four possible interactions is attractive, the other ones 
[(1, 1), (2, 1), and (2,2)] being simply HS interactions. 

In the limit of infinitely short and deep SW interactions 
(sticky limit), we discussed how a full analytical theory 
is possible. We developed a new formulation for asym- 
metric mixtures of the rational-function approximation 
(RFA), that is equivalent to the PY approximation in 
the case of symmetric SHS interactions, but differs from 
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it in the asymmetric case. Results from this theory were 
shown to be in nice agreement with MC simulations using 
SW interactions of sufficiently short width (5% of parti- 
cle size), both for the structural and the thermodynamic 
properties. 

The above agreement is rather remarkable in view of 
the rather severe approximations involved in the RFA 
analysis — that are however largely compensated by the 
possibility of a full analytical treatment — and, more im- 
portantly, by the fact that simulations deals with the 
actual Kcrn-Frenkel potential with up-down constrained 
orientations of the patches and SW attractions, while the 
RFA theory deals with the obtained asymmetric binary 
mixture and SHS interactions. We regard this agreement 
as an important indication on the correctness of the map- 
ping. 

Within the RFA approximation, all three standard 
routes to thermodynamics (compressibility, virial, and 
energy) were considered. To them we added a weighted 
average of the virial and energy routes with a weight 
fixed as to reproduce the exact third virial coefficient. 
Somewhat surprisingly, our results indicate that only the 
compressibility route fails to display a full critical behav- 
ior with a well-defined critical point. The existence of 
a critical point and a gas-liquid phase transition in our 
model (despite the fact that attractive interactions are 
partially inhibited) are supported by the NF generalized 
principle of corresponding states^ and by preliminary 
simulations results. We plan to carry out more detailed 
GEMC simulations to fully elucidate this point. 

The work presented here can foster further activities 
toward an analytical theory of the anomalous phase dia- 
gram indicated by numerical simulations of the (uncon- 
strained) Janus fluid. We are currently working on the 
extension of the present model allowing for more general 
interactions (i.e., t\i < T\\ = T22 = t%\ < 00) that can 
be continuously tuned from isotropic SHS (ry = r) to 
isotropic HS interactions (r,-j — > 00). The decrease in the 
(Baxter) critical temperatures and densities occurring in 
this case would then mimic the corresponding decrease 
in the location of the critical point upon a shrink of the 
patch width in the Kern-Frenkel models 



Appendix A: Exact low-density properties for anisotropic 
SHS mixtures 



1. Cavity function to first order in density 

To first order in density, the cavity function of an 
anisotropic mixture is 



where 



with 



(Al) 



(A2) 



V$l(?)= I dr'f ik (r')f kj (r-r>). 



(A3) 

Here, /y(r) = e~' 9 '^ ( r ) — 1 is the Mayer function. In the 
particular case of the anisotropic SHS potential consid- 
ered in this paper, 

fij(r) = /hsM + 5(r - 1) [ty0(cos0) + t Tl Q{- cosfl)] 
= ff S (r)+tr j 5( r -l)Q{cose), (A4) 

where i~ = /, , - tji, 

/hsM = -6(1 - r), /£ HS (r) = / HS (r) + t 3 -«f(r - 1). 

(A5) 



Inserting Eq. (|A4j) into Eq. (|A3|) . we get 

ygi(r) = 9(2 - r) { |(2 - r) 2 (4 + r) - (t ki + t jk )n(2 r) 



2S(r) 



1 



Hti k tjk + t kj t ki )C(v) + tr k t^C (r)} , 

where 

A(r)= { dv'5{r' - 1)6(1 - |r - r'|)6(z'), 



(A6) 



(A7) 
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£(r)= / dr'8{r'- l)5(|r - r'| -1)000, ( A ») 



C (r)= J dr' 8(r'-l)5(\r-r'\-l)Q(z')Q(z-z'). (A9) 

It can be proved that 

!7r(2-r), y/1 - r 2 /4 < costf < 1, 
A(r/2, 6), I cos<9| < ^1-^/4, 
0, -1 < cos(9 < -^/l -r-2/4, 

(A10) 
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' 2?r/r, y/1 - r 2 /4 < cos 9 < 1, 

£(r) = ^ L(r/2,0), |cos0| < y/l-r 2 /A, 
0, 



L(*) 



7T/2 



dfi sin 6>L(£,6>) 



1 < cos 6» < -^/l-r 2 /4, 

(All) 



(A22) 



2?r/r, ^1 - r 2 /4 < cos 6» < 1, 

£o(r) = <( L (r/2,e), < cos < y/l - r 2 /4, (A12) 
0, cos6><0, 



ir/2 



where 



d0 siii0L o (^0) 
i ( 7 rVl-^ 2 -2cos- 1 ^ 



(A23) 



A{£,9) = 2tt0(cos0) - tt£ - 2£sin- 1 



^ cos( 



VI - £ 2 sin* 



2. Second and third virial coefficients 



-2 tan 



.! \/sin 2 - I 2 

COS0 



1 



£ cos 6 



— h sin , 

2 v/T^Fsin0 



L o (^,0) = L(£, 6) - L(£, 7T - 6) 
2 £cos0 
£ sin 9 VT^I 2 



(A13) 



(A14) 



(A15) 



The series expansion of the compressibility factor Z in 
powers of density defines the virial coefficients: 

Z = 1 + B 2 p + B 3 p 2 + ■■■ . (A24) 

Using Eq. (|A1[) in Eq. (|4.8|) . one can identify 

B 2 = ^ 1 - 3 x i x iUj , (A25) 



In Eqs. (|A11[) and ()A12|) we have omitted terms propor- 
tional to S(r) since they will not contribute to gij{v). 
Note the symmetry relations *4(r) + A(— r) = 7r(2 — r), 
£(r) + £(-r) = 2n/r, £(r) - £(-r) = £ (r) - £ (-r). 
The orientational average 



According to Eq. (|A17|) . 



(1-3^^(1)-*^ (1) 



(i) 



(A26) 



+(*«*ifc + *«**)£"(!) + W^o(l), (A27) 



5tt 



J#*( r ) = I"'* de sm9y$ k (r) 



(A16) 



becomes 

i/giW = 6(2 - r) { 1(2 - r) 2 (4 + r) - (t ki + t jk )n(2 r) 



-tkitjk2TT 



2S(r) + - 



{^ k + t kj )A{r) 



+ t kj t ki )t{l) + tr k t kj C' (l), (A28) 

where 



— t, \ 47T /— — T l , . Z7T 

-4(l) = y-V3, .4(1) = -— 



2tt 



+ fyMr) + t- k t- 3 Co(r)} , (A17) 



where 



.A(r) = tt(2 - r) (l - y/l - r 2 /4) + A(r/2), (A18) 



£(r) = — (l - v/1-^/4) + L(r/2), (A19) 
£o(r) = — (l- -^/4) +io(r/2), (A20) 



£(l) = f, £'(1) = -H(2 7 r-V3 



_ 2n —1 2 / /-■ 

£o(l) = T , £ (1) = -3(^-2^3 



(A29) 
(A30) 
(A31) 



The second and third virial coefficients can also be 
obtained from the compressibility equation (|3.14p . To 
that end, note that 



with 



nir/2 

A{£) = I d9 $nxBA{£,i 

J sin' 1 I 



2\/l - I 2 (tt - tt£ - 1) + 2£cos" 1 £, (A21) 



where, according to Eq. (|3 . 1 3|) . 



(A32) 



(A33) 
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£g>(0) = V^72 7 r{^yg ) (l) + ^(l) 

[ drr^M+i^lV)] }. (A34) 

Inserting this into Eq. (|3.14p and making use of Eqs. 
(|A"T7)) - (|A"23"1) . one gets x^ 1 = U-2B 2 p+3B 3 p 2 + - ■ ■ , with 
£?2 and -B3 given by Eqs. (|A25|) and (|A26|i . respectively. 
Furthermore, it can be checked that the exact consistency 
cond ition (j4"X2|) is satisfied by Eqs. (jXIjl . (TMJ), fM7]) . 
and (|A28j) . The verification of these two thermodynamic 
consistency conditions represent stringent tests on the 
correctness of the results derived in this Appendix. 



3. Case tu = tii = *2i = 



yii (l) = - T + a:a* 



flV/,\ 37T 27T 

yi2 (i) = --r + *-=-, 



fll'/^x 37T 7T 

£21 ( 1 ) = -T + f 3^ 



37 (1) (1) = § -WTrt (l-f) 



(A43) 
(A44) 

(A45) 
(A46) 
(A47) 



In the preceding equations of this Appendix we have 
assumed general values for the stickiness parameters tij . 
On the other hand, significant simplifications occur in 
our constrained Janus model, where tij = tSnSj2- More 
specifically, 

yW(r) = 6(2 - r){^(2 - r) 2 (4 + r) - x 2 t [tt(2 - r) 

-tC(r)+tC (r)}}, (A35) 



|$(r)=0(2-r) -(2 - rf(4 + r) - U(r) , (A36) 



yftV) = 6(2 - r) [^(2 - r) 2 (4 + r) - x 2 7r^2 



2i r 

cos - 

7rr 2 



(A37) 



(r) = 6(2 -r) |^(2 -r) 2 (4 + r)-i ?r(2-r) 

-2a/1 -r^ + rcos" 1 r j, (A38) 



y(l)' (l ) = _^L +XlJ . 2 2t 
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S 2 =4(1 - 3xix 2 <) . 



B 3 = 10 <M - 6x10:2* 



(A48) 



(A49) 



— = —\2r\x\Xit 1 1 + - 
+0(r7 2 ). 



(A50) 

•) 

(A51) 



Appendix B: Evaluation of the coefficients \ and 

r(2) 



In order to apply Eqs. ()3.19[) and (|3.20p . it is convenient 
to rewrite Eq. (|5.12[) as 



(r) = 6(2 - r){^(2 - r) 2 (4 + r) - t [2^1^74 
}, (A39) 



— r cos — 
2 



(1) (r) = 9(2-r) ^(2-r) 2 (4 + r) -a;ia;227r^2 



* -1 r 
— cos - 
7rr 2 



m, 5?r (It 

y^ = 12 -^(1-3 



• r )7T 



(i) = ^-*hr-V5 , 
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4tt 



(A40) 



(A41) 



(A42) 



2tt 



L( S ) = Q( S )-[I-A( S )], 



where we have introduced the matrix Q as 

Qij{s) = e s s 2 G ij (s). 
Thus, Eqs. (j3~19)) and (pT2"0"|) are equivalent to 
Q i:j {s) = l + s + 0(s 2 ). 



(Bl) 

(B2) 
(B3) 



Expanding Aij(s) in powers of s and inserting the result 
into Eq. (|B1[) . one gets 



_Lr(0) _ I _ V A {0) 

k 



1 - 1 



(B4) 
(B5) 
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where 



4° = {-iTpx, 



L 



(0) 



-(i) 



(2) 



(n + 3)! (n + 2)! (n + 1)! 



(B6) 

Equations (|B4|) and ()B5j) constitute a linear set of equa- 
tions which allow us to express the coefficients and 

L\j in terms of the coefficients {L^ }. The result is 
given by Eqs. (|5.15[) and (|5.16|) . 

(2) 

It now remains the determination of L).' . This is done 



by application of Eq. (|4.17[) . i.e., the ratio first term to 
second term in the expansion of e s Gij(s) for large s must 
be exactly equal to tjj . This is the only point where the 
stickiness parameters of the mixture appear explicitly. 
The large-s behavior from Eq. (|5.12[) is 



(2) 



2^e s G i3 {s) = L i3 
where 



J#> + (L< a > • d) 



s-'+Ois- 2 ), 
(B7) 



[-Lf-Liy+L^ 



(1) , r(2) 



pxi I L 



(2) 



1-r, 



Comparison of Eq. (|4~T6l> with Eq. fB7| yields Eq. (f5TT8]) 

and 



{1 ~ >'■•.! _ r(l) 



-(2) 



fe=l 



ho- 



rn 



(BIO) 



^y(O) + hji(0) , where /ijj(O) is 



Note that h i3 (0) = \ 

defined by Eq. (|3T3]) . 

The Fourier transform /i^ (g) of the (orientational av- 
erage) total correlation function hij(r) = #«(r) — 1 is 
related to the Laplace transform Gij(s) [see Eq. (|3.15l) ] 
by 



hij(q) = -2-irpy/XiXj 



G l3 (s) - G l3 (- S ) 



(C4) 



s—iq 



Making use of Eqs. (|5.12p - (|5.14[) , it is possible to obtain, 
after some algebra, 



Cij(q) 



4tt 



00 i X i 



q 3 C if (sinq-qcosq) + -jC,; [2qsmq 



47T 



,(1) 



-2 - 
x sing + 24 
+4irUjyij(l) 



2)c0S 9 ]+£c#> [M(q 2 -6) 



- (24- 12g 2 
sing 

q 



q cosg 



(C5) 



where the coefficients , G^\ and C^' are indepen- 
dent of q but depend on the density, the composition, 
(B8) an d the stickiness parameters. Fourier inversion yields 



Cij (r) 



(o) 



y *j *j 
+^•(1)^+^-1) 



9(1 -r) 



(C6) 



Taking into account Eq. (|4.10j) we see that Eq. (|C6j) has 
the structure 



(C7) 



Cij(r) =9ij(r)-y ij (r)- 



But this is not but the PY closure relation (|5.3p . In 
passing, we get the cavity function inside the core: 



Taking into account Eqs. and ([B8]), Eq. (|Bi"0|) be- 

comes Eq. (I5.17[) . 



Appendix C: Recovery of the pseudo-PY solution 

The aim of this Appendix is to prove that the pair cor- 
relation functions (r) obtained from the RFA method 
in Sec. EH satisfy Eqs. (JESl) and (gH). 

First, note that the pseudo-OZ relation (|5.6j) can be 
rewritten in the form 



y«(r)e(l-r) = - 



c(q) = %) ■ [\ + %) 
where I is the unit matrix and 



Cij(q) = Py/XiXj / dre lq ' r Cjj(r), 



hij(q) = Py/XiXj / dre 4q ' r /ijj(r) 



(CI) 



(C2) 



(C3) 



C 



(0) 



6(1 - r). 
(C8) 
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